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Abstract. Let {x{t),y{t))~^ be a solution of a Fuchsian system of order two with three 
singular points. The vector space of functions of the form P{t)x{t) + Q{t)y{t), where 
P, Q are real polynomials, has a natural filtration of vector spaces, according to the 
asymptotic behaviour of the functions at infinity. We describe a two-parameter class 
of Fuchsian systems, for which the corresponding vector spaces obey the Chebyshev 
property (the maximal number of isolated zeroes of each function is less than the 
dimension of the vector space). 

Up to now, only a few particular systems were known to possess such a nonoscil- 
lation property. It is remarkable that most of these systems are of the type studied 
in the present paper. We apply our results in estimating the number of limit cy- 
cles that appear after small polynomial perturbations of several quadratic or cubic 
Hamiltonian systems in the plane. 

2000 MSG scheme numbers: 34C07, 34C08, 34C05 



1. Introduction 

In many bifurcation problems the main difficulty is to estimate the number of isolated 
zeroes of certain functions of the form 

I{h)=pi{h)h{h)+p2{h)h{h), /i G S, (1) 

where pi{h) and P2{h) are polynomials, and the vector function l{h) = {Ii{h), l2{h)Y 
satisfies a two-dimensional first order Fuchsian system 

l{h) = A{h)I'{h), ' = d/dh, (2) 
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with a first degree polynomial matrix A{h). Typically, Ii{h) and hih) are complete 
Abelian integrals along the ovals S{h) within a continuous (in h) family of ovals 
contained in the level sets of a fixed real polynomial H{x, y) (called the Hamiltonian) , 
and E C M is the maximal open interval of existence of such ovals 5{h). See Table 1. 
In the present paper, our main assumptions on (2) are the following: 

(HI) A' is a constant matrix having real distinct eigenvalues. 

(H2) The equation det A{h) = has real distinct roots Hq, hi and the identity 
trace A(/i) = {det A.{h))' holds. 

(H3) \{h) is analytic in a neighborhood of ho. 

The conditions that A' is a constant matrix and det A(/i) has distinct roots imply 
that the singular points of the system 

I'(/i) = A-\h)l{h) 

(including oo) are regular, i.e. it is of Fuchs type. Further, the condition trace A(/i) = 
(det A(/i))' implies that the characteristic exponents of (2) at ho and hi are {0, 1}. 
In the formulation of our main result below, we assume for definiteness that ho < hi. 
A similar result holds if ho > hi. Clearly if ho < hi, and the function I{h) is analytic 
in a neighborhood of h = ho, then it also possesses an analytic continuation in the 
complex domain C\[/ii, oo). 

Definition 1. The real vector space of functions V is said to be Chebyshev in the 
complex domain V G C provided that every function I G V^\{0} has at most dim V — 1 
zeros inV. V is said to be Chebyshev with accuracy k inV if any function I e ^\{0} 
has at most k + dim V — 1 zeros in V. 

Definition 2. Let I{h), h G C be a function, locally analytic in a neighborhood of 
oo, and s e M. We shall write I{h) < h^ , provided that for every sector S centered 
at oo there exists a non-zero constant Cs such that \I{h)\ < Cs\hY for all sufficiently 
big \h\, he S. 

For systems (2) satisfying (HI) and (H2), the characteristic exponents at infinity 
are —A and — // where A' = 1/A and = 1/// are the eigenvalues of the constant 
matrix A'. According to (H2), \-\- ji — 2. Let us denote A* = 2 if A is integer and 
A* = max (|A — 1|,1 — |A — 1|) otherwise. 

Take s > A* and consider the real vector space of functions 

Vs = {I{h) = P{h)Ii{h) + Q{h)h{h) : P,Qe R[h], I{h) < h'} 

where I = {Ii{h),l2{h)Y is a non-trivial solution of (2), holomorphic in a neighbor- 
hood oi h = ho. As ^ {0, 1, 2}, the vector function I(/i) is uniquely determined, 
up to multiplication by a constant, and Ii{ho) — hiho) = (see Proposition 1). 
Clearly, Vs is invariant under linear transformations in (2) and affine changes of the 
argument h. The restriction s > A* is taken to guarantee that Vg is not empty. 
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Recall that ho < hi are the roots of det A{h) — 0. Our main result in section 2 is 
the following. 

Theorem 1. Assume that conditions (HI) - (H3) hold. If X ^ then Vg is a 

Chebyshev vector space with accuracy 1 + [A*] in the complex domain V = C\[hi, oo). 
If X E Z, then Vg coincides with the space of real polynomials of degree at most [s] 
which vanish at ho and hi. 
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Table 1. Examples of systems for integrals I = {Ii, I2) and 
Hamiltonian functions H which satisfy hypotheses (HI) - (H3). 

As an application of Theorem 1 let us consider a polynomial perturbation of a 
planar Hamiltonian system 



X ^ Hy + ef{x,y), 
y = -H^ + eg{x,y), 
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where £ is a small parameter, the degree of the polynomials /, g does not exceed n 
and H is some of the Hamiltonians from Table 1. Define the function 



h^I{h)=l [g{x,y)dx- f{x,y)dy], /i e E. (4) 

JH=h 

As is well known, if 1(h) ^ in E, then the number of limit cycles in (3) bifurcating for 
small e from the periodic orbits of the unperturbed Hamiltonian system is bounded by 
the number of isolated zeroes of 1(h) in S. Define the linear space V„ of integrals given 
by (4) for deg f,g < n. Denote by hi the nonzero critical value of the Hamiltonian 
and by T> the complex plane cut along the part of the real axis between hi and oo 
not containing the other critical value ho — 0. Then applying Theorem 1, we obtain 
the following results. 

Theorem 2. For each of the systems 1) — 5) in Table 1, the linear space of integrals 
Vn is Chebyshev with accuracy one inV. In particular, Vn is Chebyshev in E. 

Theorem 3. For systems 6) and 7) in Table 1, the linear space of integrals V„ is 
Chebyshev with accuracy one in V, if n < 6, and with accuracy [^^]; if n > 7. In 
particular, Vn is Chebyshev in E, if n < 6, and Chebyshev with accuracy [^^], if 
n>7 



Roughly speaking. Theorems 2 and 3 imply that, for the systems l)-7) from Table 1, 
the number of limit cycles in (3) born out of periodic orbits under small polynomial 
perturbations which are transversal to the integrable directions, is less than the di- 
mension of the linear space of these perturbations (with certain accuracy if n > 7 in 
cases 6) and 7)). Clearly, a bound obtained by establishing the Chebyshev property, 
is always the optimal one. 

Case 8) from Table 1 is non- Hamiltonian one and requires slightly different ap- 
proach. See the end of the paper for results about it. 

Let us recall that Theorem 2 in case 1) was proved earlier by Pctrov [10]. Some 
less general (or a little bit different) results concerning cases 3) - 5) can be found in 
[1], [5], [6] and [8]. 

2. The Chebyshev property 

We intend first to obtain a normal form for the matrices satisfying (HI) and (H2). 
For this purpose, we perform in (2) a linear transformation bringing A' to a diagonal 
form and then translate the critical value ho to the origin. The matrix in (2) takes 
the form 

/ 2h — hi Lohi \ 
2A ~2\ 
hi 2h — hi 



A(h) 



(5) 



\ 2/100 2/1 ) 

where hi is the nonzero critical value and is a free parameter. This is the normal 
form we will use in this section. In applications, another normal form takes place. To 
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obtain it, we apply additional linear transformation in (5) {hjh) {h,h/i^ + h) 
bringing A{h) to 

/ h 



X 

(A — fi)h 

\ X/MJ 



'2X 
h hi 
H Xn J 



(6) 



Evidently, equations (5) and (6) present three-parameter families of matrices which 
can be reduced to two-parameter ones by moving hi to 1. We note that all the 
examples in Table 1 are taken in the normal form (6), with | < A < < |. 

Before to prove Theorem 1, we need some preparation. Without any loss of 
generality we may use the normal form (5), with hi — 1. Hence, we will consider 
t—{h — ho) I {hi — ho) as the argument and will assume throughout this section that 
(2) is rewritten as a system I{t) = A{t)I'{t) for I{t) = {x{t),y{t))'^ = /2(/i))^, 
with 

/ 2t-l 



A{t) 



2X 
1 



UJ 

2X 

2t-l 



\ 



(7) 



Proposition 1. The functions x{t) = Ii{h) and y(t) = hih) satisfy equations 

t{t - l)x" = A(A - l)x, 
t{t - l)y" = - l)y. 



(8) 
(9) 



Proof. The most easy proof is a straightforward calculation which we left to the 
reader (cf. [8]). 

Proposition 2. Let A 7^ 0, 1 and x{t) be a nontrivial solution 0/(8) which is analytic 
in a neighborhood oft = {ort = 1). Then x{t) ^ fort < {respectively, fort > 1). 
In particular, if X G Z, then x{t) is a special kind of ultra-spherical polynomial and 
has all of its zeros in the interval [0, 1]. 

Proof. The assertion is well known for A integer. In this case x{t) is a kind of 
ultra-spherical (Gegenbauer) polynomial [12] of degree A if A > 2 and of degree 1 — A 
if A < —1. Although the result might be known for A not integer too, we will for 
completeness give the proof for this case. Let x{t) be analytic near t — (the other 
case is similar). Take the function 



z{t) 



1 - At 



-X 



-X. 



2-A 2-A 
Then z' — tx' — Xx and x', z' together satisfy a system 



{t^ - t)x" 
{t^ - t)z" 



{X-l){tx' - z') 
{X-l){tx' -tz'). 
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As x'{0) 7^ and ^'(0) = 0, the ratio w — z'/x' is an analytical function in a 
neighborhood oi t — satisfying the Riccati equation 

^2 _ f 

-w'{t) = w^{t) - 2tw{t) + t 

A — 1 

and i/;(0) = 0. Consider in the (i, i/;) -plane the zero isocline given by the hyperbola 
up' — 2tw + t = 0. It goes through the origin and has a vertical asymptote at that 
point. It is easy to conclude that for t < 0, the graphic of w is placed inside the left 
branch of the hyperbola and either w{t) > or w{t) < for all t < 0, depending on 
whether ^'(0) is negative or positive. Therefore x'{t) and z'{t) do not change signs 
for i < 0. As x(0) = 0, the assertion follows. □ 

Proposition 3. Let A < 1 and x(t) be a nontrivial solution of (8) which is analytic 
in a neighborhood of t — 0. If X ^ Z,, then x{t) has at most 1 + [A*] zeros in the 
complex domain D — C\[l, oo). 

Proof. Consider the analytic continuation of x(t) in the complex domain D — 
C\[l,oo). We shall count the zeros of x{t) in D by making use of the argument 
principle. Let R he a. big enough constant and r a small enough constant. Denote 
by D the set obtained by removing the small disc {|t — 1| < r} from D fl {\t\ < R}. 
To estimate the number of the zeros of x(t) in D, we shall evaluate the increment 
Ag^Arga;(t) of the argument of the function x{t) along the boundary of D, traversed 
in a positive direction. Then, according to the argument principle, we have that the 
number of the zeros of x{t) in D equals 

Agf,ATgx{t) 

The monodromy group of the equation in (8) is reducible if and only if A € Z [4, 
Theorem 4.3.2]. Therefore, if A ^ Z, then in a neighborhood of i = 1 we have 

x{t)^miog{t-l) + ri{t) 

where ^{t), r]{t) are analytic in a neighborhood oit — 1, ^{t) is a non-trivial solution 
of (8), ^(1) = 0. Moreover, a local analysis shows that lim^^i- x{t) = ri{l) = const ^ 
0. Therefore the increase of the argument of x{t), when running the boundary of 
{|t — 1| < r}, is close to zero. Along the half line (1, oo) the imaginary part of x(t) 
equals Tr^{t) which does not vanish, by Proposition 2. Finally if \t\ is sufficiently big 
then we have 

\x{t)\ < c\t\^ a X > I, \x{t)\ < c\t\^-^ if X < I, \x{t)\ < c\thogt\ if X = I, 

where c is a non-zero constant. The increase of the argument of x{t), when running 
the boundary of {\t\ < R} is close to 27rA*. Summing up the above information, we 
obtain that the increase of the argument of x{t), when running the boundary of D, 
is at most 2n + 27rA*. We conclude that x{t) has at most 1 + [A*] zeros in D which 
completes the proof of Proposition 3. □ 
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We also need a more detailed information about the structure of the linear space 
Vg and an explicit formula for dim V^. The only interesting case is when A and /i are 
not integer. 

Proposition 4. Let s > X* and A, /i be not integer. Then 

{2s — 1, if A — /X and s — | are integer, 
[s — X] + [s — jj] + 2, otherwise. 



Proof. Without loss of generality, we can use the coordinates in which A takes a 
form (7) and l{h) = {x{t),y{t)y . To reduce the number of cases, let us assume that 
X> II (when X < /i, the analysis is similar). 

We begin our analysis with the case when s > A. Assume first that A — // is not 
integer. Then one can take any solution of (2) near infinity in the form 

where 

2cj(^-A + l)' 2{X- jj + l)' 

Since I is analytic in a neighborhood of zero, the constants a and b are both nonzero. 
Indeed, if a6 = then I defines an one-dimensional subspace in the space of all 
solutions, which is invariant under the monodromy group of (2), and hence of (8), (9). 
This is however impossible, as the latter groups are irreducible for A, /i ^ Z. 

Given s > A, then the function I{h) in the definition of Vg contains monomials 
of the form t^x, < k < K, t''y, < I < L, where K < min (s — A, s — + 1), 
L < min {s — X + 1, s — /i). Using that X + /i — 2 and X > /i, one obtains K < 
s - A + min (0, 2A - 1) = s - A. Similarly, L <s- ii + min (0, 2// - 1) = s - A + 1 if 
X — II > 1 and L < s — II otherwise. 

Among these monomials, other special combinations may be involved in Vg if 
A — /i > 1. Define the functions Zi = ty — aix, Zm = tZm-i — a-mX, m > 2, where 
ai = a and the constant am is determined so that the coefficient at in Zm is zero. 
Denote M ^ [s - ii] - K - 1. Clearly, then t^+'^Zm eVgior 1 <m < M. Moreover, 
any combination t^~^^{P(t)x + Q{t)y) which belongs to Vg is a linear combination of 
the "monomials" t^^^Zm- 

Thus, dim K = -f^ + + 2 for | A - /i| < 1 and dim Vg = K + L + M + 2 otherwise, 
which yields dim Vg = [s — X\ + [s — n] + 2 m. both cases. 

Assume now that A — // is integer but A and // are not. If A — // > 1, then one can 
take any solution of (2) near infinity in the form 

As in the previous case, this yields K = [s — A],L = [s — A+1] = K + 1 and 
M — [s — n] — K — 1 \l s ~ \ is not integer, M — s — n — K — 2\ls — \\s 




7 



integer. In the first case we obtain the same result as above, and in the second case 
dim K = + L + M + 2 = 2s - 1. 

Finally, if A = |, = |, we have respectively 



_3_ 
8w 



% i'' +(-— logt + 6) 1 «\ 1 



Clearly ii' = [s — |], L = [s— |] ifs — | is not integer and L — K otherwise. Since 
no other combinations are involved in Vs in this case, the result follows immediately. 

In the case when A > s > A*, the analysis is simpler. We use the same formulas 
for I as above. One has either (a) A > s > A — 1 and A — /i>l,or(b)A>s>;U and 
X — /J, < 1. Ifs = /x=|, then Vg is empty. In all other cases, y E Vg. In case (b), Vs 
contains no other functions. In case (a), if s > /x + 1 and A — // is not integer, then 
also Zm & Vg ioT 1 < m < [s — ij]. The same is true if A — is integer but s — | is 
not. Finally, if both X — ji and s — ^ are integer, and s > + 2, then V^, contains the 
functions Zm, 1 < m < [s — /i] — !. Clearly, in all the cases above we obtain a formula 
for dim Vg as asserted. □ 

Proof of Theorem 1. For integer X, fi the assertion is obvious since Ii and I2 are 
different ultra-spherical polynomials which have no common zeros except the simple 
ones at ho and hi. 

Assume below that A, ^ Z and let A > /x (for definiteness) . Suppose as before 
that the matrix A takes the form (7), and let I{t) — P{t)x{t) + Q{t)y{t) e Vg, where 
{x{t) = Ii{h), y{t) = hih)) is the holomorphic solution of (2) vanishing at the origin. 
When P{t) = 0, the assertion is evident. When P ^ 0, we use again the argument 
principle to count the zeros of I{t) in D = C\[l, 00). Consider in D the meromorphic 
function 

F{t) = P{t)^ + Qit). 

Below we calculate the increase of its argument when running the boundary of D. 
The local structure of the solutions of (8), (9) in a neighborhood of t = 1 implies that 
lim^^i x{t) 7^ 0, lim^^i y{t) 7^ 0. Therefore the increase of the argument of F{t), when 
running the boundary of {|i — 1| < r} is close to zero. As x{t), y{t) are real-analytic 
on (—00, 1), then along the half-line (1, 00) 



x(t) V x(t) x(t) 




y{t) " 2i\y{t)\' 



As {x{t),y{t)Y is the analytic continuation of {x{t),y{t)y along a loop contained in 
and the monodromy group of (9) is not reducible for /i ^ Z, then the solutions 
{x{t),y{t))~^ and {x{t),y{t))^ are linearly independent. This together with y{t) ^ for 
h G (1, 00) (Proposition 2) shows that the imaginary part of F{t) has at most degP 
zeros on (1,(X)). Suppose finally that \t\ is sufficiently big. As \y{t)\ > c|t|^* then 
F{t) < t^~^* . Summing up the above information, we obtain that the increase of the 
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argument of F(t), when running the boundary of D is at most 2n{l + deg P + s — A*). 
Moreover, in the exceptional case when A* = |, one has F(t) ~ ct*^2/logt for large 
\t\, which yields a stronger result: the increase of the argument of F on \t\ = R is 
strictly less than 27r(s — |). Therefore the total increase of the argument in this case 
is < 27r(l + deg P + s — |). This fact is useful only if s — ^ e N but we need it below. 
One can deduce from the preceding proof of Proposition 4 that: degP = [s — A] 
if A — /i < 1, degP = [s — fj, — 2] if A — /i > 1 and s — ^ arc both integers, and 
degP = [s — n — 1] otherwise. (If degP < 0, one takes P = 0.) On its hand, 
A* = /iifA — //<! and A* = A — 1 otherwise. Therefore, by Proposition 4, the 
difference between the number of zeros and poles in D of the meromorphic function 
F{t) = I{t)/y{t) is bounded by dimV^ — 1. By Proposition 3, this yields that I{t) 
has at most [A*] + dim Vg zeros in D. Theorem 1 is proved. □ 

3. The applications 

In this section we prove Theorems 2 and 3. Before that, let us point out that some but 
not everything included in Table 1 is an evident fact. However, since the procedure 
of deriving the related Fuchsian systems is more or less known, we are not going to 
discuss in more details how all these systems were obtained. 

Given i,j nonnegative integers, denote Iij{h) — jjjj^j^x^y^dxdy. Then 

Vn = {I{h)= ^'^Mh)}- (10) 

0<i+i<n-l 

Lemma 1. Let I = {hjh) be as in Table 1. Then for n > 3 one can express the 
function I{h) from (10) in the form I{h) = a{h)I\{h) + f3{h)l2{h) where a{h) and 
(5{h) are polynomials of degrees as follows: 

(i) deg a = [^], deg/3 = [^] in cases 1) and 2); 

(ii) deg a = [^], deg/3 = [^] in case 3); 

(iii) deg a = [^], deg/? = [^] in cases 4) and 5); 

(iv) dega = deg/3 = [^^] in cases 6) and 7). 

Moreover, the coefficients in a{h) and P{h) may take arbitrary values, except in case 
(iv). The dimension of the vector space Vn in the case (iv) equals [^^] + [^^] + 1- 

Proof. For some of the cases, the results in Lemma 1 arc ah^cady known. The 
result in case 1) was proved by Petrov in [5]. For cases 4) and 5) see Petrov [8], [6], 
respectively. The result for 2) follows from the considerations in [2] and [3]. The 
result concerning 3) is proved in [1]. Let us consider cases 6) and 7) from Table 
1. By symmetry, we have Iij{h) = Iji{h) and Iij{h) = whenever i or j is an odd 
number. To establish the relations between the integrals Iij{h), we take the equation 
H = x^ + y^ + vx'^y^ = h, i> = ±1 and multiply both sides by the one- form x''y^^^dx. 
Afterwards integrate the result along the oval H — h and apply Green's formula. One 
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obtains the relation 

(j + l)/,+2,i + (j + 3)J,,,+2 + l^U + 3)/i+2,i+2 = (j + l)hlij. 

Similarly, multiplying by x^^^y^dy and integrating, we get another relation 

{l + 3)/i+2,i + {i + l)-^ij+2 + + 3)/i+2j+2 = (i + l)hlij. 

Combining these equations we easily obtain 

u{i - j)Ii+2,j+2 = {j + l)^i+2,i + l)^i,i+2, i + j, 
v{i + 3)/i+2,i+2 = -(2z + 4)/i+2,i + (^ + (11) 
+ 5)/i+4,o = + 2)/i - l]i"j+2,o - 3/i,2 + hlifs. 
For i = j = 0, we get (noticing that /qo = — /i and /20 = -^02 = —^2) 

I22 = 1 1^/2 - 

Then, using (11) with i,j even, we easily prove the assertion in (iv) by induction. 

It remains to calculate the dimension of the vector space V„. Clearly, we have 
dimVi = 1, dimVa = 2, dimVs = 4. By (11), the only new functions in V2m+i 
(compared to V2m-i) are hmfi and, if m is even, Ira,m- Hence, the integrals hkfi, 
< k < m and l2k,2k, 1 < < m/2 form a basis in V2m+i- These integrals are 
independent, since the leading term of hkfi, A; > 2 is proportional to h^~^{vli + 2I2) 
and the leading term of l2k,2k, k > 1 is proportional to h'^Ii. The above argument 
implies that dimV„ = [^] + [^] + 1. □ 

Remark 1. In cases (i)-(iii) of Lemma 1, the result remains true even for n — 1,2, 
under the convention that a polynomial f3{h) of negative degree is taken to be zero. 
In case (iv), one has to take /3(/i) = 0, dego; = for n = 1, 2. 

Corollary 1. The dimension of the vector space Vn, n > 1, related to arbitrary 
polynomial perturbations of degree n in (3), in the cases 1) — 7) of Table 1 is as 
follows: 

n, in cases 1) and 2); 

[^], incases); 
2[^] + 1, in cases 4) and 5); 
[^] + [^] + 1> «™ cases 6) and 7); 

Proof of Theorems 2 and 3. Let us first note that |A — /x| < 1 for all cases l)-7) 
in Table 1, which yields that [A*] = 0. We put 

s — in cases 1) and 2), 
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s = I + I in case 3), 

s — [^^] in cases 4) and 5), 

s — [^^] + I for n > 3, s = 1 for n = 1, 2 in cases 6) and 7). 

It is easy to check that, with this choice of s, C Vs. For this purpose, one can 
perform the inverse transformation (Ji, I2) (/i, I2 — Ii/oj) bringing A to a normal 
form (5) and then use the formulas for the solution given in the proof of Proposition 
4. Hence, it suffices to verify in each case that deg a + A < s, deg (5 + n < s (the first 
inequafity should be strong in cases 6) and 7); see also Remark 1). Then we compare 
the dimensions of V„ and Vg (Proposition 4 and Corollary 1). One obtains that 

dimVn = dimV^ in cases l)-5), as well as in 6) and 7), provided that n < 6, 

dim\4 — dim V„ = [^^] in cases 6) and 7), if n > 7. 
Thus, the results follow from Theorem 1, taking into account that I{h) has always a 
zero at hg = 0. □ 

Some other examples. Let us consider in brief system 8) from Table 1. Instead of 
(3) and (4), we have 

x = Hy/M + ef{x,y), ^^^^ 
y = -H^/M + eg{x,y), 

where H — x~^{y'^ — 2x^ + x), M{x) — x"^, /, g are polynomials of degree at most n, 
and 

I{h)=(f M{x)[g{x,y)dx- f{x,y)dy], /ieS = (-l,0). (4') 

JH=h 

Note that in case 8), I is analytic in a neighbourhood oi — —1. Define by Vn the 
linear space of integrals (4') and let P = C \ [0, 00) . 

Theorem 4. For system 8) in Table 1, the linear space of integrals V„ has a dimension 
n+1 and is Chehyshev in V, with accuracy as follows: one for n — 2, two for n — 1, 3, 
three for n — and n — 3 for n > 4. 

Taking n = 2, we get the following result about the number of limit cycles in (3'): 

Corollary 2. For any quadratic perturbation of the reversible quadratic system (3'), 
the cyclicity of the period annulus around the center at {x,y) — (1,0) is two. 

Proof of Theorem 4. Denote Iki — Jjjj^^M{x)x^y^dxdy, Ik — Jjj^^M{x)x^''^ydx] 
thus I = (/2,/i)^. By symmetry, Iki = for Z odd. In the same way as above, we 
obtain the relations 

h,^ — 2fc+3f+3 (-^fc+2,; — Ik+i,i), k — —1, 0, . . . , Z = 0, 2, . . . , 

{k - \)hh+2 = (4 - 2k)h+i + (A) - 1)4, A; = 0, 1, 2, ... , 

and use the first of them to get the expression 

n 

= '^^^CkIk{h)., Cfc independent, (12) 

fc=0 
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and then the second one to obtain 



I{h) = Poh{h) + Pi{h)l2{h), for n = 0, 1, 2, 

I{h) = h-^[Piih)hih)+P2ih)l2ih)], forn = 3,4, (13) 

1(h) = h^-^[Pn-s{h)h{h) + Pn-2{h)l2{h)i for n > 5, 

where Pk denotes a polynomial of degree k. By (12), dimV„ — n + 1. Given n, 

we choose s = | if n = 0,1,2, s = ^ if n = 3,4, s = n — | if n > 5, and 
consider the corresponding linear space Vg. Its dimension is dim V? = 3, dim Vn = 5, 

4 4 

dim\/^_5 = 2n — 3, respectively. For each n the function I{h) in (13), multiplied by 
an appropriate power of h, belongs to the respective Vg. The result then follows from 
Theorem 1. □ 

Our last example is concerned with the Hamiltonian 

H = x^ + y^-x^- ax'^y^ -y^, a> 2, (14) 

which comes from the cubic Hamiltonian vector field having a rotational symmetry 
of order 4. In complex coordinates z = x + iy, such a field is presented by a complex 
equation z = —iz + Az'^z + Bz^, A, B E C, Re A = 0. Take a polynomial perturbation 
in (3) which is semi-even with respect to x: 

f{-x,y) = f{x,y), g{-x,y) = -g{x,y), degP,Q<n (15) 

and consider the integral (4) where E = (^^) |) ^^^^ the integration is along the oval 
S{h) C {H = h} surrounding the center at (^^j^O). As in Lemma 1, we can derive 
relations between the integrals involved in (4) and then use them to rewrite I{h) in the 
form I{h) = P{h)Ii{h) + Q{h)l2{h) where h = f^^^-^ x'^dy, h = J^^f^^ x^y^dy, and P, Q 

are polynomials with independent coefficients and degrees [^^], [^^], respectively. 
The related vector space Vn has a dimension [|]. The vector function I = (Ji, 12)^ 
satisfies a system (2) with a matrix (which is too large to fit in Table 1) 

/ 4/1-1 a -2 \ 

_ I 3 3 

4/i - 1 4/i a-U 
\ 15(a + 2) 15(a + 2) / 

Clearly, conditions (H1)-(H3) are satisfied with ho = ^, hi = Denote V — 

C\ (— 00, /ii]. Take s = then evidently Vn — Vg. Applying Theorem 1, we obtain 

Theorem 5. For any system (3) satisfying (14) and (15), the linear space of integrals 
Vn has a dimension [|]. Moreover, Vn is Chebyshev with accuracy 1 in V and it is 
Chebyshev in E. 

Theorem 5 is useful for estimating the number of limit cycles not surrounding the 
origin that are born in small semi-even polynomial perturbations of the cubic Hamil- 
tonian vector field with a rotational symmetry of order 4. 
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Appendix: Non-oscillation and Sturm type theorems 

The classical Sturm theorem can be used to find bounds for the number of the zeros 
of the solutions of linear non-autonomous differential equations on a real interval. In 
the context of the present paper a Sturm type non-oscillation theorem was recently 
proved by Petrov [8] . The proof uses of course topological arguments. It is natural to 
ask whether the results of the present paper couldn't be deduced in such a way. The 
answer turns out to be negative in general, and our main Theorem 1 is essentially 
a non-oscillation result in a complex domain. On the other hand our proofs also 
rely on topological arguments: the argument principle for real analytic functions in 
a complex domain. Therefore we may call Theorem 1 a Sturm type theorem in a 
complex domain. 

To compare these two approaches (real and complex) we give below an example 
in which a Sturm type theorem in a real domain can still be proved. We shall follow 
closely Petrov [8]. As in the introduction, P and Q are the polynomials from the 
definition of Vg and it is assumed (for definiteness) that ho < hi. 

Theorem. Assume that conditions (H1)-(H3) hold and2X ^ Z. Then any nontrivial 
function in Vg has at most degP + dcgQ + 1 zeros in the interval (— cxD,/to)- -^^ 
particular, if \X — /i] < 1, then Vg is a Chebyshev vector space in (— oo, ho). 

Proof. As in section 2, it is sufficient to consider (2) as a system for 1(h) — 
{x{t),y{t))~^ , with A taken in a normal form (7). For k a nonnegative integer, denote 

(^2k+i = ik + \){k + X - 1), uj2k+2 = {k + (i){k + n - 1), 

V UJ2k+2 J V '^1^ ii + k-l J 

Following [8] we introduce the operator 

By Proposition 1, we have LI = VLqI. Next, we prove that under hypotheses (Hl)- 
(H2), the operator L satisfies also the following identities: L(t'^I) = t^Q,k^ — i'^^^Rfcl, 
A; G N. Indeed, taking into account the form of the matrix in (7), and denoting for 
short 5 — det A = t{t — 1)/XijL, we obtain 

Mt'^i) = t{t - i){t''i)" 

= t(t - l)[t^r + 2kt^-^r + k{k - l)t^-H] 

= [t'^L + 2kXfit''-^SA-^ + k{k - l)(t^ - t''-^)]l 

= t''[Qo + 2kXn{6A-^y + k{k - 1)]I + t''-^[2kXfi{SA-^){0) - k{k - 1)]I 
— t^Vt]^ — t^ ^Rj^I- 
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Assume that 2A is not integer. Then it is easy to verify that the constants cuj are all 
different. This implies that there exist scalar functions of the form 

y,{t) = 0{t^-^)x{t) + [t^ + 0(t^-i)]y(t), 
satisfying the equations 

LXkit) = UJ2k+lXk{t), Lyk{t) = U!2k+2yk{t), 

where 0(t'^^^) denotes different polynomials of degree k — 1. To verify this, we ask 
for a Ifc(t) = {xk{t),yk{t))~^ in the form 

k 

i=o 

Bfc the unity matrix, to be determined for j < k. As the operator L commutes 
with the constant matrices, we have 

k-l 

Llfc = flkt'^i + — Bj+iRj+ijt-'I = flk^k 

3=0 

and the matrices Bj, j — k — 1, k — 2, . . . ,0 are determined recursively from the 
equations 

Hjflj — OfcBj = Bj+iRj+i 

which is possible because Uj are all different. Therefore, there is a basis in the space 
Vg consisting of the eigenfunctions of the operator L. Taking into account that, by 
Proposition 2, and do not vanish for t < 0, we apply Petrov's elimination 
technique [8] to prove that any function in Vg has at most deg P + deg Q + 1 isolated 
zeros in (— cxd,0). Especially, in the case when |A — /i| < 1, this means that Vs is 
Chebyshcv in (— oo,0). □ 

Note that the above proof works only on the open intervals having ho as an end- 
point and where det A is positive (because the Sturm theorem applies in a backward 
direction here). Also note that when |A — /i| > 1, the above estimate, although it 
concerns the interval (— oo, /lo) only, is weaker than the estimate for the whole V 
obtained in Theorem 1. 
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